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ABSTRACT 

Recent multi-waveband observations of Seyfert nuclei and QSOs established significant deviations 
in the spectral shape of the big blue bump from a blackbody one; Soft X-ray excess has a spectral 
index a (F^ cx v^") of 1.6 and hard X-ray tail with a of ^ 0.7. We construct a disk-corona model 
which accounts for such broad-band spectral properties. We study emission spectrum emerging from a 
vertical disk-corona structure composed of two-temperature plasma by solving hydrostatic equilibrium 
and radiative transfer self-consistcntly. A fraction / of viscous heating due to mass accretion is assumed 
to be dissipated in a corona with a Thomson optical depth of Tc, where advective cooling is also included, 
and a remaining fraction, 1 — /, dissipates within a main body of the disk. Our model can nicely reproduce 
the soft X-ray excess with a power-law shape and the hard tail extending to ~ 50 keV. The different 
spectral slopes (a ~ 1.5 below 2keV and ~ 0.5 above) are the results of different emission mechanisms 
and different sites; the former slope is due to unsaturated Comptonization from the innermost zone 
and the latter is due to a combination of the Comptonization, bremsstrahlung and a reflection of the 
coronal radiation at the disk-corona boundary from the inner to surrounding zone (< 300 Schwarzschild 
radii). The emergent optical spectrum is redder {a ~ 0.3) than that of the standard disk {a ~ —0.3), 
being consistent with observations, due to the different efficiencies of spectral distortion of disk emission 
at different radii. Further, we find that the cut-off frequency of the hard X-ray coronal electron 
temperature) and broad-band spectral shape are insensitive to the black-hole mass, while the peak 

frequency of the big blue bump is sensitive to the mass as the peak frequency cx A/bh^^^^- 



Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: nuclei 
— radiation mechanisms: miscellaneous 

1. INTRODUCTION 

Active Galactic Nuclei (AGNs) exhibit a lot of spectral components over a wide waveband; optical/ultraviolet (UV) 
bump, power-law component in hard X-rays, excess component in soft X-ray band (so called soft excess), warm absorber 
feature around 1 keV in some objects, and reflection/Fe Ka fluorescence line around a few to a few tens of keV (Mushotzky, 
Done & Pounds 1993; Koratkar & Blaes 1999 for reviews). Those spectral components are commonly thought to be 
powered by gas accretion onto a massive black hole. The most well-known disk-accretion model is called as the standard 
disk model (Shakura & Sunyaev 1973). According to the standard disk model, the spectrum at each radius of the disk is 
assumed to be a blackbody radiation with a local effective temperature, T^g. This simple picture was supported by their 
luminosity and by rough agreement of optical/UV spectral energy distribution (SED) between observations and models 
(e.g.. Shields 1978; Malkan & Sargent 1982). 

The standard model, however, has limitations; (i) if UV turn-over around 2000 A is indeed an indication of the 
temperature at the innermost radius (e.g., Malkan 1983; Sun & Malkan 1989) such disks are too cool to produce enough 
soft X-ray photons, (ii) Inversely, if the soft excess component at 0.1 to 1 keV is due to radiation from the innermost 
region, as is often interpreted, such a disk only produces by 1.5-2 orders of magnitude less optical/UV flux than what is 
observed, (iii) When optical and soft X-rays spectra are simultaneously fitted with the disk spectrum, the luminosity of 
the disks (L) often exceeds the Eddington luminosity (iEdd)- (iv) Observed optical spectra of QSOs are typically redder 
{a ^ 0.3, where L^, oc i/^"; Francis et al. 1991) than those of the simplest standard accretion disks (a ^ —1/3). In other 
words, a successful model spectrum of AGNs should deviate from that of the standard accretion disk, (v) Hard power-law 
X-ray cannot be reproduced. 

Then, a number of authors have tried to distort the disk spectrum toward the high energy regime so that the disk can 
emit substantial soft X-ray and optical/UV radiation simultaneously. One promising idea is the Comptonization within 
the disk in the vertical direction (Czerny & Elvis 1987; Wandel & Petrosian 1988; Laor & Netzer 1989; Ross, Fabian & 
Mineshige 1992). The most accurate treatment of Comptonization in the framework of the standard model was made 
by Shimura & Takahara (1993, hereafter ST93) who solved radiative transfer and vertical structure simultaneously and 
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presented emergent spectra integrated over radii (see also Shimura & Takahara 1995, hereafter ST95). The effect of 
Comptonization is more prominent at higher accretion rates (e.g., Ross et al. 1992; ST93; ST95). 

However, there still remain discrepancies between models and observations, (i) Although Comptonization tends to 
increase a, the Far-UV (FUV) spectra of these accretion disks exhibits a ~ 1 at best when the luminosity of the disk 
approaches the Eddington limit (Ross et al. 1992; ST95). On the other hand, the observations of distant quasars showed 
steeper FUV spectra (a ~ 1.8 - 2.2, Zheng et al. 1997). (ii) The observed small spectral indices in soft X-ray (a ~ 1.4 
- 1.6; e.g., Walter & Fink 1993; Laor et al. 1997) are not achieved by any disk models (e.g., Nandra et al. 1995 for 
Mrk 841; Laor et al. 1997 for a sample of low-redshift QSOs), because the higher energy tail of those disk models is the 
superposition of Wien laws, thus exhibiting exponential roll-over, (iii) These accretion disks still cannot reproduce the 
hard X-rays. Then, hard X-ray emission should be treated as additional components in these models. 

With these problems kept in mind, we, in this present study, aim to produce the overall SED simultaneously by disk- 
corona models. The spectrum to be reproduced with these models is a composite one obtained from several independent 
observations (see Zheng et al. 1997; Laor et al. 1997). It is composed of 'typical' spectral indices over broad-bands from 
Near-IR to hard X-ray, optical-to-X-ray flux ratio, and the energy cut-off of the hard power-law component (see §3.1 
for more detailed description). We can, for the first time, reproduce such broad-band spectra. In §2 we review the 
assumptions and numerical method used in the calculation. Numerical results and comparisons between the models and 
the composite spectrum are presented in §3. The final section is devoted to discussion and summary. 

2. BASIC ASSUMPTIONS AND EQUATIONS 

The numerical code used in this study is basically the same as that of ST93 except for some modifications. We 
assume for the configuration of the system that the accretion disk main body is sandwiched between coronal layers in the 
vertical direction (e.g., Haardt & Maraschi 1991), and that the whole system is geometrically thin (i.e., plane-parallel slab 
geometry) . We treat the disk-corona system consisting of fully ionized hydrogen, thermal plasma around a Schwarzschild 
black hole of mass Mbh- Gas evaporation from the disk to corona and -condensation from corona to disk (e.g., Meyer & 
Meyer-Hofmeister 1994) are not included here for simplicity; the disk and coronal layers are interacting only via radiation 
and pressure (e.g., Nakamura & Osaki 1993; Zycki, CoUin-Souffrin & Czerny 1995). Although magnetic fields may also 
affect the disk-corona structures, they are not included in this calculations. 

The equation of hydrostatic equilibrium in the vertical direction is given by 

where R and z are the radial and the vertical coordinates, respectively; and Pgas and are the gas pressure and the 
number density of electrons, respectively. For the radiation field, we adopt the diffusion approximation. Then, the 
radiative fiux at some frequency v is given by 

where e^, and A^, are the radiation-energy density per unit frequency per unit volume and the mean free path of a photon 

with a frequency z/, respectively. We set A,y = (ctt -I- as_u)~^ we do not include the bound-free nor free-bound 
transitions in the calculations. Here, ctt and trff^^ are the cross sections of Thomson scattering and bremsstrahlung 
absorption, respectively. The Thomson optical depth of the disk-corona measured from the mid-plane r is related to the 
height from the mid-plane z, as 

dr = Ne (TT dz. (3) 

The total Thomson optical depth from the mid-plane (through the boundary between the disk-corona) to the surface of 
the corona, tq, is a free parameter. The equation of state is that Pgas = N^k-QiTp + Tg), where Tp and Tg are temperature 
of protons and electrons, respectively. 

Th{^ dissipated energy per unit surface area of the disk-corona system Qq in Newtonian approximation is written as 
(Shakura & Sunyaev 1973) 

Q+ = ^Mn^ [i - ^mjR)j , (4) 

where M and Ok [= (GA/ehZ-R"^)^^^] are the total (above and below the mid-plane) accretion rate and the Keplerian 
angular velocity, respectively. The innermost radius of the disk-corona is assumed to be 3iig, where Rg [= 2GMbh/c^ = 
3 X lO^^cm(AfBH/lO^Af0)) = 0.01 It-day (Mbh/IO^^M©)] is the Schwarzschild radius. For the case of a non-rotating black 
hole imder non-relativistic treatment, M — 12LEdd/c^ [= 2.6 M0yr^^(AfBH/lO*MQ)] corresponds to the accretion disk 
shining at LEdd [= 1-3 x lO'^^erg sec-^(MBH/lO^M0)]. 

A constant fraction / of mass accretion is assumed to be dissipated in the corona with a Thomson optical depth of Tc, 
where advective energy transport of protons is also included in addition to radiative cooling of electrons (see Figure 1). 
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A remaining fraction, 1 — /, dissipates within the disk layer. Advection in the disk layer (i.e., slim disk model that will 
be discussed later) is not included because of small radial velocity there. The advective cooling rate in the corona per 
unit surface area, Q~^^, is taken from the expression of optically-thin advection dominated accretion flow (ADAF) in 
one-temperature case (i.e., Tp = Te; see, Kato, Fukue & Mineshige 1998, p272): 



167r^ J acnipTc 



O- - — — — ^ ^ f'Sl 



Here, etc is the viscosity parameter which controls the efficiency of the advection in the corona. We perform numerical 
calculations under the condition that the advective cooling should be less than the dissipated energy in the coronal layer; 

Qa"dv < / Qt- 

The heating rate in the disk and corona, and , respectively, and advective cooling rate in the corona per unit 
volume, q~^^, are assumed to be proportional to matter density at each site. As a result, the fraction of advective cooling 
in the dissipated energy is 

^^dv _ Qa-dv _ o.2i_i_ I __ 1 ( _^ ) li-W::^) . (6) 



(it fQo OicTc \L-Edd/C- 

The energy balance in each layer is as follows: 

Disk : Protons 
Electrons 
Corona : Protons 
Electrons 




= Aie (7) 

Aie = ^rad (8) 

it = ^adv + (9) 

Aie = «-d' (10) 



with 



= {^- f)Qt——^ Qt=fQo^—^ 9adv=Qadv-— . (11) 

'0 'c 'c 'c 

where Aic is the energy exchange rate due to Coulomb collisions taken from Guilbcrt and Stepney (1985); and q^^^ is 
the radiative cooling rate. We consider bremsstrahlung and Comptonization for emission mechanisms. Then, radiative 
cooling rate ^ is described as 

9rad = j dl'{-cNeaf[^^e^+jf[^^+jcomp,u), (12) 

where js^i, and icomp,!/ are the bremsstrahlung emissivity (Rybicki fc Lightman 1979) and the net rate of energy transfer 
from electrons to photons via Comptonization per unit volume per unit frequency, respectively. The latter is described 
by the Kompaneets equation (Rybicki & Lightman 1979; see Hua & Titarchuk 1995 for comparisons of the analytical 
treatment of Comptonization with Monte Carlo simulations). Thus, the radiative transfer equation for F^, is written as 

-g^ = -cNeas,^^^ + + jComp,v (13) 

For the expression of the free- free absorption and emission [cra^v and js,v, respectively), we take the Gaunt factor to be 
unity, for simplicity (see Rybicki & Lightman 1979). 

Following ST93, we take ^ [=log(To — r) ] as an independent variable for the vertical coordinate. To calculate the spectra 
integrated over the whole disk (§3.1), we divide the disk-corona from 3i?g to 300i?g into 20 consentric rings so that each 
ring radiates approximately the same luminosity (cf. Ross et al. 1992; ST95). In total, input free parameters required for 
the calculations are Mbh, M, f, Tc, tq, and ac- The number of these parameters is similar to that of relevant observed 
parameters which we aim to explain simultaneously; e.g., Lx, cvox, ct\jv, ctopt, ctROSAT, ciasca, etc. Then, outputs are 
A^e(C)) Te{£,), Ti{(,), Pgas(C)i ^(Oj and e^(^). Spectrum of the whole disk-corona system is obtained by summing up the 
emergent spectra of all the rings. 

The emergent spectrum at each ring does not so strongly depend on tq as long as tq ^ 1, but is sensitive to Tc- We, 
thus, take tq as a constant over all rings, for simplicity. Currently, we do not have a good theory to predict the radial 
dependence of Tc and / that can, in principle, be determined by physics of evaporation/condensation. Thus, Tq and / are 
also simply assumed to be constant over all rings. The effects of changing Tc and / will be discussed later. 

General relativistic effect is not included. Since no Doppler broadening is considered, the total spectrum represents the 
case of a face-on disk. Convection/conductive energy transport is not included. Shakura, Sunyaev & Zilitinkevich (1978) 
have shown that convection transports no more than 30 % of the vertical energy flux. 
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3. SPECTRAL ENERGY DISTRIBUTION OF AGNS 

3.1. Observed Composite Spectrum 

The spectrum to be reproduced with our models is a composite one (dotted lines in Figure 2a), which is a useful probe 
to see whether a model for the SED of AGNs can work or not. The vertical normalization is determined so as to give 
rise to a representative optical luminosity among low-redshift quasars from the Bright Quasar Survey adopting Hq of 50 
km s^^ Mpc^^ and go of 0.5 with an assumption of isotropic emission (sec Laor et al. 1997). There arc, however, two 
major problems that one should keep in mind when dealing with the composite spectrum (Koratkar & Blaes 1999; Laor 
1999). (i) We observe objects with different redshifts in different wavebands; spectral index in FUV is obtained from 
distant QSOs (Zheng et al. 1997), while that of soft X-ray is mainly from nearby objects (e.g., Walter & Fink 1993; Laor 
et al. 1997). Also, sample used in each wavebands contains objects that are not necessarily the same objects, (ii) Soft 
excess might be an instrumental/calibration problem; BeppoSAX observations did not detect soft X-ray excess in some 
objects while excess was seen for the same objects with other telescopes (e.g., Matt 1999). Although the results are still 
controversial, we assume that soft excess really exists in most AGNs throughout this paper. In what follows, we try to 
reproduce broad-band spectra similar to the composite one, as the first step. 

3.2. Model Spectrum Integrated over Radii 

We first show the most successful case with Mbh = 3 x IO^A'/q and M = 0.5 LEdd/c^ that gives rise to liuninosity of 
about 5 % of the Eddington limit for a non-rotating black hole under non-relativistic treatment. Thick line in figure 2a 
shows an example of the resultant broad-band spectra. Since the innermost ring from 3.0 Rg to 4.0 Rg does not satisfy 
the restriction that (J^civ/^c^ — ^-^ ^'-'^ ^^"^ parameter sets described in Fig. 2a, wc, hereafter, plot integrated spectra from 
4.0 Rg to 300 Rg. In other words, all dissipated energy is assumed to be carried out by advection in 3.0 < R < A.ORg. To 
account for the radiation at the surfaces of coronae above and beyond the disk, the resultant spectrum is multiplied by 
two. A significant fraction / = 0.6 of mass accretion occurs at the corona. In this case, advective cooling in the corona is 
comparable with the radiative cooling at iZ ~ 4.9i?g; q^^/qt ^ 0.5 [cq. (6)]. A spectrum of the standard disk with the 
same Mbh and M is also depicted for comparison (dashed line). 

Presence of multiple spectral components is the most noteworthy feature of the present model. This is because dif- 
ferent radiative mechanisms play roles in different wavebands in Fig. 2a; thermal radiation of the disk in optical/UV, 
unsaturated Comptonization in FUV/soft X-ray, and a combination of the power- law component due to Comptonization, 
bremsstrahlung, and a reflection in hard X-ray. Note that the underlying radiative processes in soft -hard X-rays are 
distinct from those of the traditional explanation, in which UV soft X-ray component is due to blackbody whereas hard 
power-law component is due to Comptonization of the soft photons. In our model, hard X-ray spectrum looks a power-law 
with a ~ 0.5-1.0 due to the combination of multiple radiative mechanisms, in contrast. 

Figure 2b shows contributions to the total spectrum (4.0-Rg < R < 300i?g; thick line) due to individual rings. The 
outermost ring (104i?g < R < 300i?g) contribute as much as one third of the total spectrum at a few tens of keV. 
Usual models for X-ray emission of AGNs assume that only inner region radiates X-ray (i.e., one-zone); our model differs 
from the traditional model in terms of radial dependence of the X-ray spctrum, as well. Like optically-thin ADAF (e.g., 
Narayan, Yi & Mahadevan 1995; Manmoto, Mineshige & Kusunose, 1997), X-ray emission arises from wide spatial range 
{R ^ 300i?g), in contrast with a usual belief. Interestingly, contribution from the outer rings to the total X-ray spectrum 
decreases towards lower-energy X-ray band. Provided that inner rings are more time-variable than outer rings, resultant 
X-ray spectrum will get softer when the luminosity increases, as is actually observed (e.g.. Done et al. 1995; Leighly et al. 
1996). The most powerful test of our model will be gravitational microlensing (Yonehara et al. 1998) which provides 
information as to the size of emission region as a function of wavelength on AU (~ Rg for IO^Mbh) scales. 

The eS'ect of advective cooling in the coronal layer to the emergent spectrum is demonstrated in Figure 2c. Thick 
dashed line is the resultant spectrum without advective cooling, which corresponds to a two-temperature treatment of the 
study by Shimura, Mineshige, & Takahara (1995), but with a constant Tc- Lower dashed curves represents the spectra of 
each ring in the case without advective cooling. Thick solid curve has the same meaning as in Fig. 2a and 2b. 

In our models above (Fig. 2a-2c), we assumed constant / and Tc over all rings. This will be a reasonable assumption 
since Janiuk, Zycki & Czerny (2000), for example, estimated / as a function of radius under some assumptions, finding 
that / is a slowly increasing variable with radius, and that 0.2 < / < 1.0 for various parameter sets at 3 to 300 Rg. The 
coronal gas is probably originated in evaporated gas cumulated in the flow (Meyer & Meyer-Hofmeister 1994; Liu, Meyer 
& Meyer-Hofmeister 1997). In other words, we have assumed that a disk corona has already formed until ~ 300 Rg. If, 
however, evaporation of disk material is still active at i? ^ 300i?g, Tc and / are both likely to increase with a decreasing 
radius. To what extent the resultant spectrum changes in such a case is demonstarated in figure 2d. Spectra for the 
outermost ring with f = Tc = 0.6 (solid curve) and f = Tc = 0.1 (dashed curve) are shown. In the latter case, hard X-ray 
emission is not strong as the former case, in which bremsstrahlung emission at the outermost ring is the main origin of 
hardening in the total X-ray spectrum (thick solid curve). Then, we are lead to the conclusion that the corona must have 
developed until 300i?g so as to have large values of / and Tc at the outermost ring, thereby reproducing the observed 
spectrum. 



3.3. Red Optical Spectrum 
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Optical spectral index of QSOs has been recognized as a discrepancy between disk models and observations; Francis 
et al. (1991) reported that the bright QSOs, in which stellar contamination to the optical fluxes is relatively negligible, 
typically show spectral index in optical, a at 1500-5000 A, of 0.32 [L^ oc while the classical standard disk model 

predict that of —0.33 (i.e., oc u^^^). 

In order to emphasize the optical spectra, a part of Fig. 2a is enlarged and shown in figure 3 with additional spectra 
of disk models integrated over 3 < i? < lOOORg. Thick solid and dotted lines have the same meanings as in Fig. 2a, and 
they both show similar spectral index, a ~ 0.3. The shaded area inndicates the relevant waveband, 1500-5000 A. 

The reason why the current model has a better agreement with the observations in terms of the optical spectral index 
is as follows. The long-dashed line is a spectrum of the standard disk (i.e., eq. 4), but without the inner boundary term 
(the term within the parenthesis on the right-hand-side of eq. 4, which is sometimes omitted when roughly estimating 
the spectral index of the disk model). It shows a ~ —0.3. Spectra of the standard disks with the inner boundary term 
are drawn by the short-dashed lines. The upper one represents the disk with Mbh = 3 x lO^M© and M = 0.5LEdd/c^) 
while the lower one is for Mbh = lO^M© and M = Z/Edd/c^ showing somewhat redder spectrum {a ~ —0.1) due to the 
inner-boundary term. 

The emergent local spectrum changes its shape in a sense that the spectrum is distorted and shifted towards shorter 
wavelength when the electron scattering in the disk is taken into account (Czcrny & Elvis 1987). Such a spectral shift is 
more efficient in inner region than in outer region. As the result, disk models with electron scattering finally show redder 
optical spectra, a <^ +0.1 in lower solid curve and a ~ +0.3 in upper solid curve. The difference of spectral indices in 
the two solid curves, which are obtained by radial integrations of ST93 model (i.e., no corona), arc due mainly to the 
difference of their black-hole masses. The disk with larger Mbh gets cooler at fixed M/iEdd and R/Rg, as is expected 
for the standard accretion disk (eq. 14 described later). The optical band becomes to measure the flux at inner part of 
the disk, where it gets close to the edge of the integral interval (e.g., 3 < i? < 1000i?g). Then, the emergent spectrum 
is somewhat curved in a sense that it is redder than the spectrum from outer region (e.g., figure 16 in Koratkar & Blaes 
1999). Note that the redder spectra are achieved even within the framework of original ST93 and ST95 models, as well 
as our disk-corona model. If one wants to discuss in more details about optical spectra, heavy metals, bound-bound, free- 
bound, bound- free transitions and precise estimation of the contribution from broad line region (e.g., Balmer continuum) 
must be included. 

3.4. Spectra for Various Mbh and M 

Figure 4 shows the black-hole mass dependence of the spectra with a fixed ratio of M to iEdd (i-c, fixed M/Mbh)- 
The peak frequency of the big blue bump varies with Mbh in a similar fashion to that of the standard accretion disk in 
which 

Teff(4i?g) oc Mbh-'/'M^/* oc MBH-'/^M/LEdd)'/*. (14) 

On the other hand, the cut-off frequency of the hard X-ray (i.e., coronal electron temperature) and the spectral slopes in 
X-rays are rather insensitive to the black-hole mass. 

Accretion-rate dependence of the emergent spectra with a fixed blackhole mass (Mbh = 3 x IO^Mq) is shown in Figure 
5. The spectral slope at 0.1 - 2.0 keV stays almost constant, a = 1.46, 1.58, 1.62 for M/(LEdd/c^) = 0.1, 0.5 and 0.7, 
respectively. 

3.5. Vertical Structure 

We finally present the vertical structure of the disk/corona. Number density of electrons (dotted line) and temperature 
of protons (crosses) and electrons (solid line) at i? ~ 4.9i?g are shown in Figure 6. The left side of the figure where ^ [= 
log(ro — t)] =3 corresponds to the mid-plane (i.e., tq = 1000), and the boundary between the disk and corona is located 
at ^ —0.2 (i.e., Tc = 0.6). Parameters used here are the same as those of Fig. 2a. It turns out that the height of the 
boundary measured from the mid-plane is 0.03 Rg, and that of the surface of the corona (at ^ of —2) is 0.3 Rg. Then, 
the disk-corona system is indeed geometrically thin; the height of the corona from the mid-plane ^ O.IR. 

4. DISCUSSION AND CONCLUSION 

4.1. Comparison with previous models 

The most widely accepted interpretation, to date, of soft and hard X-ray emission mechanisms is that soft excess is due 
to unscattered photons propagating from the disk-corona boundary to the coronal surface, while hard power-law emission 
is attributed to unsaturated Comptonization, i.e., photons that are Compton up-scattered during the propagation within 
the corona. Monte-Carlo simulations of propagating photons in the framework of two zone treatment (disk and coronal 
layers) indeed reproduce bump-like feature in soft X-ray and power-law component due to Comptonization in hard X-ray 
(e.g., Haardt & Maraschi 1991; Nakamura & Osaki 1993). However, the bump-like feature due to unscattered photons 
still keeps its spectral shape as it was at the disk-corona boundary even after it goes away from the coronal surface: i.e., 
it should look like blackbody radiation. Then, higher energy tail of the soft excess is predicted to be as steep as the Wien 
law, in contradiction with observed small spectral indices in soft X-ray (a ^ 1.4 - 1.6). 



6 



Broad-Band SEDs of AGNs from an Accretion Disk with Advective Coronal Flow 



Thermal bremsstrahlung emission was proposed for hard X-ray power-law component by Schlosman, Shaham & Shaviv 
(1984), in which the emission comes from the disk-corona transition layer with an temperature gradient, where thermal 
conduction from the overlying corona is balanced with radiative cooling. They show that such a layer emits hard X-rays 
efficiently in the case of an accretion disk model proposed by Sakimoto & Coroniti (1981), in which viscous torque is 
proportional to gas pressure, but that it does not work in standard accretion disks where the torque is in proportion 
to the total pressure unless viscosity parameter in the disk is less than 10^"^. We should note here that non-thermal, 
power-law electron energy distribution models, similarly to our results, produce two spectral power-laws for soft and hard 
X-ray emissions (e.g., Zdziarski & Lightman 1985), although ours is a thermal model. 

Shimura et al. (1995) proposed a model in which FUV-hard X-ray is attributed to unsaturated Comptonization in a 
corona above a disk. Their claim is that their model, with Tc being a strongly increasing function with increasing radius, 
reproduces a single power-law in X-ray bands, not a broken power-law (i.e., soft excess and hard power-law). Then, the 
corona could be patchy so that the disk is covered over only in part with the corona. In our model, however, bremsstrahlung 
emission from the corona and reflection of the coronal emission at the disk-corona boundary also contribute to hard X-ray 
as well as the unsaturated Comptonization. As a result, we have different spectral slopes at soft and hard X-ray bands. 
It largely depends on our assumption of constant / and t^. In other words, we need / and to be almost constant or 
weak functions of radius in order to expect sufficient bremsstrahlung emission from the corona at outer region and then in 
order to reproduce the observed spectrum within the current framework. Theoretical works to check the assumption are 
needed as next steps. Without advective cooling in the corona, we also have rather straight spectrum from FUV to hard 
X-ray (see Fig. 2c). The effect of advection seems to bent the total spectrum so as that the spectrum has two power-laws 
in X-ray bands. 

In the extremely high accretion rate, the disk spectra may be able to reproduce the observed steep spectral decline 
at 1000A-400A, since Comptonization is more efficient at higher accretion rates and since an increase of accretion rate 
changes the disk dynamics from the standard disk to optically-thick ADAF. Optically-thick ADAF, which is so called the 
slim disk model, radiates with higher temperature than the standard disk (Szuszkiewicz, Malkan & Abramowicz 1996; 
Mineshige et al. 2000). Narrow-Line Seyfert 1 galaxies (NLSls) have peculiar spectral and temporal features that are not 
seen in normal Seyfert nuclei and QSOs (e.g.. Boiler, Brandt, & Fink 1996; Brandt, Mathur, & Elvis 1997; Leighly 1999a, 
1999b; Grupe et al. 1998, 1999). These features are often attributed to small black- hole mass and to high accretion rate 
(i.e., L ^ L-Edd] e.g.. Pounds, Done, & Osborne 1996; Hayashida et al. 1998; Mineshige et al. 2000). Then, such extreme 
accretion rate is not relevant in this study, where we are going to construct a disk model for normal Seyferts/QSOs that 
perhaps have moderate accretion rate [L ~ (0.01 - 0.1) of LEdd; e.g., Wandel 1999]. To account for systematically large 
oiRosAT in NLSls, we would need to reduce Compton-y parameter in the corona. 

The configuration of cold (~ 10^ K) and hot (~ 100 keV) regions that radiate optical/UV bump and X-rays, respectively, 
is an unsolved problem. Optically-thin ADAF is a possible energy source of X-ray emissions. However, optically-thin 
ADAF itself is basically faint (e.g., Mahadevan 1997) compared with observed luminosity of AGNs. Then, the luminosity 
of the power-law X-ray emission indicates that cold material is located adjacent to such hot flow as a source of seed 
photons: for example, corona-like flows above and below a disk (e.g., Liang & Price 1977; Haardt & Maraschi 1991). 
Separation in the vertical direction has been also suggested from the stability against thermal and secular instabilities in 
inner, radiation-pressure dominated region of the disk (lonson & Kupcrus 1984: Nakamura & Osaki 1993; Mineshige & 
Kusunose 1993), the correlated variabilities between the iron-line flux and the X-ray continuum (e.g., Nandra et al. 1999; 
see also, however, Lee et al. 1999), and need for a cold disk as a 'mirror' for reflection-component/broad-fluorescence line 
(e.g., Tanaka et al. 1995; Mushotzky et al. 1995; Nandra et al. 1997). Another potential problem is that those layers 
could be inhomogeneous and highly time dependent; the corona can be patchy (Shimura, Mineshige, & Takahara 1995; Di 
Matteo 1998; Kawaguchi et al. 2000; Machida et al. 2000), or the cold gas may exist as numerous blobs within the corona 
instead of the slab geometry (Guilbert & Rees 1988; Lightman & White 1988; Sivron & Tsuruta 1993; Collin-Souffrin 
et al. 1996; Krolik 1998; Rozariska 1999). Nevertheless, spatially and temporally averaged treatment of the cold and hot 
regions will still be useful for study of time-averaged spectrum. Appropriate treatment of mass evaporation/condensation 
process is another future issue. In other words, / and Tc may be variable in terms of R, Mbh, -M, etc. For instance, see 
Janiuk et al. (2000) for / as a function of R with various M and viscosity parameter, and Zycki et al. (1995) for Tc v.s. 
M. 

4.2. Expected Time Variability 

Time variability provides strong constraints on spectral models. Then, it will be worthwhile to briefly comment on 
the temporal behaviour, although we are now working on a model for the steady accretion disk-corona with an aim to 
reproduce the observed, time-averaged spectrum. Followings are qualitative arguments expected from our current model. 
More details are out of the scope of this paper, and will be remained as future works. 

For the optical-to-UV time-lag observed in NGC 7469 (Wanders et al. 1997; Collier et al. 1998) and marginally detected 
in NGC 4151 (Peterson et al. 1998), our model can explain the lag as Collier et al. (1998) demonstrated using the standard 
accretion disk model. That is because optical/UV emitter in our model is, like other spectral models, an accretion disk 
that has a similar dependence of (surface) temperature upon radius to the standard accretion disk. It should be noted 
here that such a time-lag has not yet been established in all AGNs (Edelson et al. 2000) 

Since the FUV, soft X-ray and a part of hard X-ray emission come from the same region in our model, they are expected 
to vary without serious time-lag (larger than an order of days), except for difference of escaping (diffusion) time at different 
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energies. Chiang et al. (2000) reported the inter-band lags in NGC 5548 from ~ 240 ksec simultaneous observations; 
0.14-0.18 keV flux leads 0.78 keV flux by ^ 10 ksec, and 5.4 keV flux by ~ 30 ksec. As discussed in section 6 in their 
paper, these lags are consistent with a coronal size of 10 Rg for a 10® M© black-hole. Incidentally, the blackhole mass 
of NGC 5548 is estimated to be about 10*Mq through the reverberation mapping method (Wandel, Peterson & Malkan 
1999; Kaspi et al. 2000). However, the ^ 4-day X-ray lag behind UV near the peak flux levels and almost simultaneous 
variations near the minimum fluxes detected in NGC 7469 (Nandra et al. 1998) is not solved by our model easily, as well 
as other models, to date, fail to do. 

Finally, our model qualitatively explains following two issues about spectral variability in X-rays on short timescales 
(less than months) that are reported by numerous X-ray observations (§3.2). i) Soft X-ray flux is more variable than hard 
X-ray flux, ii) Hard X-ray becomes softer when it gets brighter. 

4.3. Conclusions 

We study emission spectrum emerging from the vertical disk-corona structure composed of two-temperature plasma 
by solving hydrostatic equilibrium and radiative transfer self-consistently. The key question is what physical condition 
exhibits a soft X-ray excess with a spectral index a [F^, cx u~°') of 1.6 and a hard X-ray tail with a of 0.7 at the same 
time. In our model, a fraction / of viscous heating is assumed to be dissipated in a corona with a Thomson optical depth 
of Tc, where advective cooling is also included, and a remaining fraction, 1 — /, dissipates within a main body of the 
disk. Our model can nicely reproduce the observed composite spectrum of AGNs, which shows soft X-ray excess with a 
of about 1.5 and hard tail extending to ~ 50 keV with a different slope {a ^ 0.5). Our model should be checked with 
individual objects in the future, though. 

The broken power laws {a ^ 1.5 below 2keV and ~ 0.5 above) are the results of different emission mechanisms: 
unsaturated Comptonization in soft X-rays and a combination of the Comptonization, bremsstrahlung, and a reflection of 
the coronal radiation at the disk-corona boundary in hard X-rays. Previous models, where soft X-ray excess is attributed 
to blackbody or saturated Comptonization of the disk blackbody, tended to deal with limited wavebands separately, while 
we tried to fit the broad-band SED simultaneously. That is the reason why we propose here emission mechanisms, that 
are different from previous models, in soft and hard X-ray bands. Also, our model differs from traditional models in 
terms of X-ray emitting region; hard X-ray emission causes from spatially spread region up to 300 Schwarzschild radii. 
The emergent optical spectrum is redder {a ^ 0.3) than that of the standard disk (a ~ —0.3), being consistent with 
observations, due to the different efficiencies of spectral distortion and shift of disk emission at different radii. The cut-off 
frequency of the hard X-ray (reflecting the coronal electron temperature) and X-ray spectral slopes are insensitive to the 
black-hole mass, while the peak frequency of the big blue bump is sensitive to the mass. 
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Fig. 1. — A schematic view of the energy balance in the present model. 
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Fig. 2a. — Resultant spectrum from the disk-corona structure integrated over 4.0 to 300i?g (thick line). Parameters used in this 
model are listed in the figure. With those parameter sets, advective cooling in the corona is comparable with the radiative cooling at 
R ~ 4.9Rg; q~^^/qt ~ 0.5. Dashed line indicates the integrated spectra of the standard disk with the same Mbh and M. Spectral 
indices of the observed composite spectra (dotted lines) arc: a of 1.4 at NIR (A > l^m), 0.3 at optical (l^m > A > 2500A), 1.0 at 
UV (2500A > A > lOOOA), 1.8 at FUV (lOOOA > A), 1.6 at soft X-ray (0.2-2.0 keV), 0.7 at hard X-ray (> 2.0keV), 1.5 between 
optical and X-ray (2500 A-2keV). The energy cut-off for the hard X-ray power-law is assumed to be 100 keV. 
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Fig. 2b. — Contributions to the total spectrum due to individual rings (lower curves). Thick line denotes the total spectrum 
(4.0J?g < R< 300J?g) as the same as in Fig. 2a. 




Fig. 2c. — The effect of advective cooling in the corona layer to the emergent spectrum. Thick dashed line is the resultant 
spectrum without advective cooling. Lower dashed curves represents contributions to the total spectrum due to each ring. Thick 
solid curve has the same meaning as in Fig. 2a and 2b. 
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Fig. 2d. — Spectra for the outermost ring with f = Tc = 0.6 (soUd curve) and f = Tc = 0.1 (dashed curve). It is shown that hard 
X-ray emission in the latter case is not strong as the former case, where brcmsstrahlung emission at the outermost ring is the main 
origin of hardening in the total X-ray spectrum (thick solid curve; see Fig. 2b). 
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Fig. 3. — Spectra of various models with optical spectral indices a {L^ oc z/~™). The shaded area inndicates the relevant waveband, 
1500-5000 A. Thick solid and dotted lines have the same meanings as in Fig. 2a. Spectra of the standard disks [Shakura & Sunyaev 
1973 (SS73), i.e., eq. 4] are shown as the short-dashed lines as in Fig. 2a. A long-dashed one is calculated through eq. 4 without 
inner boundary term (the last term in eq. 4). Finally, thin solid curves are obtained by radial integrations of ST93 model (i.e., no 
corona). Lower three calculations are made for Mbh = 10* M0 and M = LEdd/c^ at 3 < R < WOORg, while upper two curves are 
for Mbh = 3 IO^Mq and M = 0.5 iBdd/c^ at the same radii. 
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Fig. 4. — Black-hole mass dependence of the emergent spectra with a fixed ratio of M to LEdd- Parameters used here are the 

same as those of Fig. 2a and 2b. The cut-off frequency of the hard X-ray (reflecting the coronal electron temperature) and X-ray 
spectral slopes are insensitive to the black-hole mass, while the peak frequency of the big blue bump is sensitive to the mass. 
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Fig. 5. — Accretion-rate dependence of the emergent spectra with a fixed black-hole mass (Mbh = 3 X 10^ Mq). 
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Fig. 6. — Vertical structure of the matter density (dotted line), proton temperature (crosses), and electron temperature (solid 
line). The mid-plane of the disk is located at ^ of 3, and the boundary between the disk and corona is at ^ —0.2. It turns out 
that the height of the boundary measured from the mid-plane is 0.03 Rg, and that of the surface of the corona (at ^ of —2) is 0.3 
Rg. Note that k-BT/m^c^ = 1 corresponds to T ~ 5 10'' K (~ 500 keV). 



